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ABSTRACT 

A  class  of  polychotomous  quantal  response  models  is  defined,  including 
the  dichotomous  response  model.  As  an  application  of  the  Missing  Information 
Principle,  it  is  shown  that  the  Fisher  information  for  the  common  parameters 
in  the  more  elaborate  model  is  greater.  For  the  logistic  curve  and  several 
design  schemes,  the  above  result  is  supported  in  a  numerical  study.  For 
estimating  the  common  parameters  of  the  dichotomous-response  and  trichotomous- 
response  models  and  the  percentiles  of  the  common  tolerance  curve,  an 
extensive  simulation  study  shows  that  in  general  the  latter  model  gives  more 
accurate  estimates  than  the  former  one.  The  situations  in  which  such  gains 
are  substantial  are  identified. 
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SIGNIFICANCE  AND  EXPLANATION 


A  dichotomous  quanta 1-response  model  is  commonly  used  to  relate  the 
response  probability  to  the  underlying  stimulus  variable.  The  percentiles  of 
the  quantal  response  curve  are  often  used  as  measures  of  the  quality  of  the 
test  item.  Sometimes  it  is  practically  feasible  to  classify  the  responses 
into  three  (or  more)  ordered  categories,  say,  animal  alive,  moribund,  or 
dead.  A  basic  issue  is  whether  any  (or  how  much)  gain  in  statistical 
efficiency  can  be  realised  from  this  more  elaborate  classification  of  the 
response.  We  give  a  positive  answer  both  theoretically  and  empirically,  and 
identify  the  situations  in  which  such  gains  are  substantial. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC ,  and  not  with  the  authors  of  this  report. 


COMPARISON  OP  THE  DICHOTOMOUS  AND  POLYCHOTOMOUS  QUANTAL  RESPONSE  MODELS 
Siu-Keung  T»e  and  C.  P.  Jeff  Hu 

1.  Introduction 

In  quanta 1  response  model,  the  observed  response  is  impossible  to  measure 
quantitatively  but  instead  is  a  nominal  variable  which  usually  takes  on  two  possible 
outcomes,  response  or  non-response.  This  modsl  is  characterized  by  assuming  a  tolerance 
distribution  P  for  every  subject  in  the  population.  Thus,  if  stress  level  x  is 
applied  to  a  test  subject,  then  the  response  y  is  a  random  variable  taking  value  1 
with  probability  P(x)  and  O  with  probability  1  -  F(x).  nils  model  is  common  to  many 
areas  of  research.  In  engineering  the  stress  levels  may  ba  the  pounds  force  applied  to 
equipments  or  the  temperature  increases  in  a  system.  Then  the  number  of  brittle  failures 
or  unsuccessful  performances  can  ba  recorded.  In  material  testing,  glass  containers  are 
submitted  to  a  drop  height  test,  the  response  is  "broken*  or  "not  broken*.  In  drug 
testing,  if  the  houseflies  are  exposed  to  a  particular  insecticide,  the  number  of 
houseflies  that  are  either  dead  or  alive  can  be  observed. 

Usually,  the  tolerance  distribution  assusws  the  form  of  a  probit  or  logistic  model. 
Some  other  models  include  the  angular  distribution  model,  the  rectangular  model,  etc. 
Finney  (1978)  compared  the  various  models  and  concluded  that  they  are  indistinguishable 
between  response  rates  0.05  and  0.95.  Thus,  in  our  numerical  and  simulation  study,  we  use 
the  logistic  function  for  the  tolerance  curve  but  the  theoretical  results  hold  more 
generally. 

However,  in  the  examples  mentioned  above,  there  may  be  more  refined  classifications 
of  the  responses.  Por  example,  a  response  can  further  be  classified  according  to  the 
varying  degrees  of  damage  or  moribundity  besides  the  two  definite  responses.  Aitchison 
and  Silvey  (1957),  Ashford  (1959)  examined  the  problem  based  on  the  maximum  likelihood 
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2.  Comparison  of  Informations  in  Largs  Samples 
2.1.  Polychotomous  Quantal  Response  Model 

Assume  that  k  groups  consisting  of  N1>N2> . • . ,Nk  subjects  are  tested  with  stimulus 
levels  x1,x2>. .. respectively.  Out  of  the  subjects  being  tested  at  level  x^, 

the  possible  ordered  outcomes  of  the  test  are  ,02, .. . ,0g.  In  the  housefly  example, 
for  s  ”  3,  0|  ■  dead.  02  “  moribund,  03  »  alive.  For  a  >  3,  the  0^'s  represent 
decreasing  degrees  of  moribundity  as  i  increases.  For  consistency,  the  first  outcome 
represents  dead  and  the  last  alive.  For  simplicity  and  clarity,  we  use  this  housefly 
example  to  describe  the  model. 

For  level  x^,  let  r^j, . . . ,ria  be  the  numbers  of  subjects  with  outcomes 
Of,...,Os  respectively,  i.e.  •  r^  +  ri2  +  •••  +  ris.  The  observed  proportions  of 
outcomes  are  p^j  -  i  m  1,2,. ...a,  i  “  1,2, ...,k.  Assume  the  outcome  Oj  has  a 

tolerance  distribution  function  of  the  form  F(x.0,«j)  where  F  satisfies: 

i)  The  first  and  second  partial  derivatives  of  F  exist  and  are  not  identically 
zero, 

ii)  F(x,0,a)  is  a  strictly  monotone  function  of  a. 

Note  that  0  is  restricted  to  be  constant  for  all  the  tolerance  distributions. 
Therefore,  with  assumption  (ii),  the  tolerance  curves  would  not  cross  each  other  so  that 
the  probability  of  each  outcome  Oj  is  well-defined.  Gurland,  Lee  and  Dahm  (1961) 
pointed  out  this  fact.  For  instance,  the  logit  and  probit  models  satisfy  (i)  and  (ii). 

tfithout  loss  of  generality,  we  assume  a,,  <  a2  <  •••  <  Therefore 

-  P{xi  >  tolerance  level  for  Oj) 

»  F(xif0,aj)  j  -  1,2,. ...s  -  1  . 

Hence  P{ outcome  is  °1>  "  Pi1 

Ptoutcome  is  °j>  ■  PD  - 
p{ out come  is  0^}  “  1  -  Pj^.^  • 


J  -  2 


s  -  1 


The  piru>t«ri  in  these  models  ara  $  .a,  ,a2»  •  •  •  *®s_i  •  Our  main  interest  is  to  estimate 
the  percentiles  of  the  tolerance  function. 

Define  the  EDp  or  LDp  to  be  the  level  Xp  such  that  on  the  avera9e  100p 
percent  of  the  test  subjects  will  die  at  this  dose  level.  Then,  we  have 

p  -  FUp.B.a,)  • 

Therefore,  the  estimation  of  EDp  depends  only  on  the  parameters  6,a^.  The  log- 
likelihood  function  of  this  model  takes  the  form 


log  La  -  ^  fri1lo9  Pi1  +  ri2log(Pi2  ~  pit>  +  ***  + 


<Ni  -  ri1  .  ri.s-l110?'1  -  Pi,s-1)}  * 

T 

Denote  9  ■  (P,a1 ,a^, . . . •aB_1 )  •  The  Fisher  information  matrix  is  *  (J^jlsH, 
1  <  i,j  <  s  where 

2 

3  logL  3 logL  3  logL 

Ji:)(s)  -  E(  ae  ^  39  ^  ’  ”E^  39  39*^  ' 

i  1  i  3 


Under  some  regularity  conditions  on  F(x,0,a),  the  matrix  Js  is  well-defined  and 
is  positive  definite.  Furthermore,  the  asymptotic  distribuiton  of  the  HLE  of  6  is 
N(e,JB’),  i.e.  the  asymptotic  covariance  matrix  is  the  inverse  of  the  Fisher  information 
matrix.  Therefore,  the  Fisher  information  matrix  is  a  measure  of  the  size  of  the 
asymptotic  confidence  ellipsoid. 


One  can  partition  Jg  into 
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where  J . 


is  an  m  x  m  matrix,  Ji 


for  any  m  *  1,2, 


an  (s  -  m)  x  (s  -  m)  matrix.  In  particular. 


3 logL  3 logL  3 logL  3 logL 
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-32logL  -321o$L  -32logLs 

303q2  iWflj 

-3 2  logL  -32logL  -32logL 

. . a  •  •  •  ^ 
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n>«n  the  Fisher  information  matrix  corresponding  to  (B,a,)  is  J 


(2)  _ 
11 


t(2),(2)_1_(2) 


'12  u22 


J2 1 


2.2.  Dichotomous  Quanta 1  Response  Model 

If  only  two  outcomes  are  recorded,  either  dead  or  alive,  then  this  model  is 
equivalent  to  combining  the  outcomes  Oj.Oj, ...,08  into  one  class,  namely,  the  total 
number  of  subjects  that  are  not  dead. 

With  the  notations  of  $1,  the  proportions  are  pi1  ”  dead,  and 

1  -  plt  -  (st  -  r^J/t^  not  dead.  Similarly, 

K(pii>  *  P^,  ”  P{xA  >  tolerance  level  for  0,} 

-  F(x1,6,«1)  . 

The  Fisher  information  matrix  for  (8,0,)  under  this  model  is 

f  -32log|L2  -32logL21 


J2 


3B2 

32logL2 

383a, 


3B3a, 

32logL2 


(4) 


2.3.  Comparison  of  the  Fisher  Information  Matrices 

In  a  (s  +  D-reaponse  model,  the  numbers  of  subjects  with  outcomes  0,  ,02, . . .  ,0S+, 
at  level  x1  are  *'it»ri2' •  •  • 'ri, s+1  respectively.  Let  Sj  ■  (Ni,ri, , . . .  ,ri> s+, ) , 
i  ■  1,2,...,k  be  the  vectors  of  responses  for  the  polychotomous  response  model.  If  the 
outcoaws  0-#0B+,  are  collapsed  into  one  class,  then  the  number  of  response  subjects  in 
this  combined  category  is  the  sum  of  the  corresponding  rij,s.  Equivalently,  Os  and 
°s+,  can  viewed  as  missing  data  with  only  the  total  Os  +  Os+,  available.  Therefore, 
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th«  information  available  in  tha  lattar  a-reaponaa  modal  can  ba  aussaarlxed  by  the  vectors 
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The  proof  is  given  in  the  Appendix. 


By  a  similar  argument,  we  can  show  that  more  information  is  gained  by  having  a  more 
refined  classification  of  other  ordered  responses.  Therefore,  although  there  are 
additional  parameters  to  be  estimated  in  the  augmented  model,  there  is  still  an 
information  gain  as  a  result  of  refining  the  classification  of  the  response. 

In  particular,  the  next  result  compares  the  Fisher  information  matrix  associated  with 
the  parameters  (B,a1 )  of  any  polychotomous  response  model  with  that  of  the  corresponding 
dichotoanus  response  model . 

Theorem  2i  Consider  the  Fisher  information  matrix  Js  of  an  s-response  model  and 
partition  Jg  into 

where  jj2*  is  the  2*2  upper-left  sub-matrix.  Then  J^2*  “  JI1'J22*  J21*  *  J2  for 

any  s  >  3. 

The  proof  is  given  in  the  appendix. 

This  theorem  states  that  the  Fisher  information  matrix  associated  with  the  parameters 
(0,aj)  in  the  polychotomous  model  is  always  greater  than  or  equal  to  the  Fisher 
information  matrix  for  the  correaponding  dichotomous  response  model.  From  this  result,  it 
is  concluded  that  more  information  is  gained  in  using  the  polychotostous  response  model  for 
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•stinting  the  EDp,  or  any  smooth  function  of  (8,0^ )  over  the  dichotomous  response 
model  in  large  samples. 

2.4.  Evaluation  of  the  Fisher  Information  Matrices 

In  this  section,  ve  perform  a  small  sample  comparison  of  the  Pisher  information 
ntrices.  In  particular,  we  assume  the  underlying  tolerance  distribution  function  is 
logistic  and  s  “  3,  i.e.  01  ”  dead,  02  “  moribund,  O3  “  alive.  Therefore  Pi1  * 
probability  of  dead  •  1/[1  +  exp(-8x^  -  o^)],  Pi2  "  probability  of  dead  or  moribund 
-  1/(1  +  exp(-8xt  -  a2)>,  a2  >  a,. 

Three  different  design  patterns  are  chosen.  Scheme  I  chooses  the  10,  40,  60,  90 
percentiles  of  the  standard  logistic  curve  (8  *  1,  ^  ■  0]  as  design  points.  Scheme  II 
chooses  the  5,  40,  60,  95  percentiles  while  the  S,  45,  55,  95  percentiles  are  used  in 
Scheme  III.  Note  that  these  levels  are  symmetric  and  the  probability  of  moribund  would  be 
affected  solely  by  a2 .  Therefore,  different  values  of  a2  will  change  the  relationship 
between  the  two  tolerance  curves.  Various  numbers  of  observations  are  allocated  to  each 
percentile.  The  design  configurations  vary  from  the  balanced  to  the  unbalanced,  either 
skewed  to  the  left  or  right. 

The  criterion  for  comparing  the  two  measurement  procedures  is  the  volumes  of  the 
corresponding  confidence  ellipsoids  of  the  parameter  estisiates  for  a  specified  confidence 
coefficient.  The  volume  of  the  confidence  ellipsoid  depends  on  the  determinant  of  the 
corresponding  asymptotic  covariance  matrix.  Therefore,  such  a  comparison  reduces  to  the 
comparison  of  the  generalised  variances  of  the  asy^totic  multivariate  normal  distribution 
involved.  Thus,  the  ratio  of  the  determinants  of  the  corresponding  asymptotic  covariance 
■wtrices,  or  equivalently,  the  inverses  of  the  Fisher  information  matrices  based  on 
different  sndels  gives  a  numerical  measure  of  the  asymptotic  relative  efficiency.  In 


Table  1,  we  computed  the  ratio  of  the  determinants  of  the  Fisher  information  matrices 
corresponding  to  (8,<>j)  under  the  two  models.  Ihe  formulae  are  given  by  (1),  (2),  (3) 
(4)  in  Sections  2.1  and  2.2. 


Table  1:  Ratio  of  the  Determinants  of  the  Fisher  Information  Matrices 


for  the  3 -Response  and  2-Response  Model  ( 8  “  t,  o1  *  0) 


Design  Scheme 

<*2  ■  0.1 

dj  -  0.5 

BB 

d2  *  1.5 

' 

<&2  -  2*0 

QQ 

d2  -  4.0 

(15,5,5,15)  I 

1.04 

1.18 

1.29 

1.32 

1.28 

1.15 

1.06 

II 

1.04 

1.22 

1.41 

1.50 

1.48 

1.28 

1.12 

III 

1.05 

1.23 

1.43 

1.53 

1.50 

1.29 

1.12 

(10,10,10,10)  I 

1.04 

1.19 

1.36 

1.44 

1.43 

1.27  !  1.12 

II 

1.04 

1.24 

1.49 

1.70 

1.79 

1.57 

..27 

III 

1.05 

1.25 

1.53 

1.77 

1.87 

1.62 

*  29 

(15,15,5,5)  I 

1.04 

1.21 

1.39 

1.48 

1.47 

1.29 

1.13 

II 

1.05 

1.29 

1.61 

1.85 

1.93 

1.65 

1.30 

III 

1.06 

1.32 

1.67 

1.94 

2.01 

1.67 

1.31 

(5,5,15.15)  I 

1.04 

1.18 

1.33 

1.43 

1.44 

1.29 

1.14 

II 

1.03 

1.18 

1.39 

1.58 

1.69 

1.56 

1.28 

III 

1.03 

1.17 

1.38 

1.58 

1.70 

1.58 

1.29 

(5,15,15,5)  I 

1.04 

1.19 

1.37 

1.50 

1.54 

1.39 

1.19 

II 

1.04 

1.23 

1.49 

1.76 

1.96 

1.90 

1.49 

III 

1.05 

1.25 

1.57 

1.92 

2.18 

2.11 

1.60 

Note:  (15,5,5,15)  I  means  that  15,5,5,15  observations  are  respectively  assigned  to  the 

four  levels  of  design  scheme  I. 


The  two  tolerance  curves  are  related  in  the  following  manners: 
i)  If  a2  is  close  to  dj,  then  the  two  curves  are  close  to  each  other,  i.e.  the 

probability  P^2  “  pii  of  being  moribund  is  small  and  most  of  the  subjects  observed 
are  either  dead  or  alive. 
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ii)  Similarly,  If  a2  la  too  large,  the  curve  corresponding  to  01  +  02  ia  cloae 
to  1,  l.e.  the  probability  1  -  P^2  of  being  alive  is  email  and  the  subjecta 
obaerved  are  mostly  dead  or  moribund. 

In  either  caae,  one  of  the  claaaes  02,0j  has  few  subjects.  Therefore,  it  is  not 
worth  spending  additional  labor  in  refining  the  response  02  *  Oj  into  02  and  Oj. 

Thus,  the  gain  in  using  the  trichotomous  response  model  in  these  cases  would  not  be  very 
significant.  As  we  can  see  from  the  table,  the  ratio  increases  as  a2  increases  to  2.0 
and  then  decreases.  In  fact,  when  a2  is  sufficiently  large,  the  corresponding  ratio  of 
the  determinants  of  the  information  matrices  is  close  to  1.  These  facts  can  be  seen  by 
observing  that  the  likelihood  equations  are  the  same  by  letting  a2  ♦  a1  or  a2  ♦  “. 

The  percentage  gains  shorn  in  the  table  are  at  least  5%  with  most  of  them  more  than 
ISt,  some  up  to  100%.  This  gives  an  asymptotic  justification  for  choosing  the  3-response 
model  to  estimate  the  parameters.  We  have  also  computed  the  ratio  of  the  determinants  of 
the  Fisher  information  matrices  for  other  design  schemes.  They  exhibit  the  same  pattern, 
and  are  therefore  omitted. 

In  the  above  comparison  we  use  the  determinant  of  the  Fisher  information  matrix  as  a 
convenient  summary  criterion.  A  practical  question  is  to  what  extent  doer  this  criterion 
reflect  the  performance  of  estimates  of  different  parameters  in  the  finite  sample 
situations?  It  can  only  be  answered  by  Monte  Carlo  simulation.  Such  results,  reported  in 
the  next  section,  support  the  general  conclusion  of  this  section  that  the  use  of  the  3- 
response  model  gives  rise  to  improvement  in  efficiency  of  estimates. 
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3.  A  Simulation  Study 


Tha  experimenters  ara  of tan  intaraatad  in  the  aatiaation  of  tha  percentixes  of  tha 
tolaranca  distribution  aaaociatad  with  death.  In  the  housefly  exaatple,  we  may  be 
interested  in  estimating  the  dose  level  such  that  100p%  of  tha  houseflies  are 

killed.  In  this  section,  w»  will  present  an  extensive  simulation  study  to  compare  the 
trichotomous  response  and  dichotomous  response  models  for  estimating  0,  a^,  ED10,  ED30, 
ED50,  ED70,  ED90  under  different  design  schemes .  The  Fisher  information  in  {2  is  an 
asymptotic  measure  while  Table  2  presents  simulation  results  for  sample  size  40.  The 
latter  results  provide  a  more  realistic  indication  of  the  performance  of  these  estimates 
in  small  samples  under  the  two  models. 

Again,  we  assume  the  underlying  tolerance  distribution  F(x,0,a)  is  logistic. 
Therefore,  the  100pt*1  percentile  of  the  distribution  function  corresponding  to  death  is 
given  by 

Xp  »  -ta,  ♦  logt  1  -  p)/p]/0  .  (5) 

A  A 

The  MLE  of  xp  is  then  obtained  by  substituting  0^,6,  the  MLE  of  0^,0  into  (5).  The 
MLE  a1 ,0  of  a,B  are  found  by  using  the  iterative  Newton-Raphson  method. 

The  design  schemes  are  the  same  as  those  considered  in  Table  1.  For  each  of  the 
design  scheme  considered,  we  generated  1000  random  samples.  The  method  of  maximum 
likelihood  is  used  in  estimating  the  parameters.  The  results  are  summarized  in  Table  2. 
The  MSE  is  computed.  The  true  values  of  the  parameters  are  ■  0.0,  0  ■  1.0.  By  (5), 
we  have  ED10  -  -2.1972,  ED30  -  -0.8473,  ED50  •  0.0,  ED70  -  0.8473  and  ED90  »  2.1972. 
Therefore,  if  <*^,0  are  the  MLE  of  a^,0,  then  the  MLE  of  EDp  is  given  by 
ED_  -  [-a,  +  ED  ]/0.  Note  that  the  MLE  a,,  ED50  of  a,  and  ED50  may  not  take  the 

PIP  1  I 

same  value  although  o1  »  ED50  »  0. 
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Table  2:  Comparison  of  MSE  of  the  KLE  of  Parameters  a^,S  and  the  Percentiles 

In  Standard  Logistic  Model  with  lOOt  Simulation  Samples  (a^  -  0,  (5  *  1  ) 


Design 

Scheme 


3-Response  Model 


Parameter 


2-Response 

Model 


(15,5,5,15)  I 


0.26  0.26  0.24  0.23  0.25  0.26  0.27 
0.33  0.35  0.32  0.31  0.31  0.34  0.34 
0.33  0.33  0.31  0.31  0.31  0.33  0.34 


0.15  0.15  0.11  0.10  0.11  0.15  0.18 


0.28  0.24  0.19  0.14  0.13  0.20  0.25 

0.38  0.36  0.29  0.20  0.19  0.37  0.53 


0.65  0.62  0.52  0.51  0.53  0.64  0.72 
0.68  0.63  0.53  0.51  0.51  0.61  0.69 
0.73  0.65  0.55  0.52  0.51  0.63  0.72 


0.26  0.25  0.22  0.22  0.23  0.26  0.28 
0.30  0.29  0.26  0.26  0.26  0.29  0.31 
0.31  0.29  0.26  0.26  0.27  0.30  0.31 


0.20  0.20  0.19  0.18  0.19  0.20  0.20 

0.24  0.24  0.22  0.23  0.23  0.24  0.24 


0.24  0.24  0.22  0.23  0.24  0.24  0.24 


0.28  0.28  0.27  0.25  0.26  0.26  0.27 
0.30  0.30  0.29  0.29  0.29  0.30  0.30 
0.31  0.31  0.29  0.30  0.30  0.31  0.31 


0.70  0.69  0.63  0.58  0.59  0.65  0.70 
0.67  0.66  0.62  0.59  0.59  0.63  0.66 
0.71  0.68  0.63  0.62  0.61  0.67  0.72 


Sunnary  of  Tabla  2: 


1.  For  naarly  all  the  design  schenes  considered,  the  estimates  of  the  parameters  (8, a,) 
from  the  trichotomous  response  model  have  smaller  MSB  than  those  from  the  dichotomous 
response  model.  The  reduction  is  about  50%  in  the  best  case.  Therefore,  even  for 
small  sample  sices,  the  results  of  42  still  hold.  This  gives  an  empirical  support  of 
those  results. 

2.  Similar  to  the  performance  of  the  Fisher  information  matrix  discussed  in  |2,  the  MSB 

of  the  parameter  estimates  8  and  a^  are  minimized  at  ■  1.5  or  2.  But,  the 
MSB  increases  when  a2  is  too  large  or  small.  Therefore,  the  trichotomous  response 
model  is  useful  only  if  the  two  curves  corresponding  to  and  01  +  o2  are  not  too 

close  to  each  other,  or  if  the  curve  corresponding  to  01  +  02  is  not  too  close  to 

1.  See  the  explanation  after  Table  1. 

3.  The  simulation  study  shows  that  there  is  a  significant  improvement  in  estimating  8 
by  using  the  trichotomous  response  model  while  the  gain  is  not  so  significant  in 
estimating  a.,.  This  can  be  explained  by  observing  that  the  slope  parameter  8  is 
common  to  the  two  curves  in  the  3-response  model  and  can  be  estimated  efficiently  by 
’borrowing  strength*  from  data  in  all  categories. 

4.  For  the  estimation  of  the  extreme  percentiles,  ED10  or  ED90,  there  is  large 
reduction  in  MSB  by  using  the  3-response  model  instead  of  the  2-response  model.  On 
the  other  hand,  the  reduction  is  not  so  remarkable  for  the  median  ED50.  a  heuristic 
justification  is  given  as  follows. 

Since  xp  «  Hog  -  o^l/B  . 


where 


var(x  ) 
P 


varfa.j) 

cov(a,,8> 


cov(a1 ,8) 
var ( 8 ) 


3x 
— E 
3a, 


3x 

_J2 

36 


3x  , 

— E  _  — 

3a1  8 
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VA.  . 

0.5 

Therefore,  when  p  -  0.5,  “  0,  8  ■  1  implies  — *  0  and 

var(x0  5>  •  var(a^)  . 

As  pointed  out  in  3)  the  estimation  of  a^  and  equivalently  of  xo.s  does  not 
benefit  much  froa  the  second  tolerance  curve  in  the  3-response  model . 

5.  When  a2  is  too  large  or  saall,  i.e  a2  *  0. 1  or  4,  the  MSE  froa  the  trichotoaous 
response  model  nay  produce  slightly  larger  KSE  than  those  froa  the  dichotoaous 
response  model.  However,  the  differences  in  all  cases  are  less  than  5%.  nils  nay  be 
due  to  the  sampling  fluctuation  in  the  Monte  Carlo  simulation. 

6.  For  the  estimation  of  low  (high)  percentiles,  the  designs  that  are  skew  to  the  right 
(left),  i.e.  allocating  more  subjects  at  the  low  (high)  dose  levels,  perform  better 
than  the  others.  The  MSE  associated  with  the  extreme  percentiles  ED 10  and  ED90, 
are  higher  than  those  of  030,  ED50,  ED70.  This  illustrates  the  fact  that  the 
estimates  of  these  extreme  percentiles  are  very  unstable.  In  particular,  the  V-ahaped 
design  (i.e.  more  observations  on  the  extresw)  gives  the  smallest  MSE  for  these 
extreme  percentiles  while  those  associated  with  the  K-pb+.'A  design  (i.e.  fewer 
observations  on  the  extremes)  are  the  largest.  However,  the  balanced  design  gives 
overall  good  performance  in  estimating  the  percentiles. 
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4.  Conclusion 


Generally,  our  sain  cones rn  is  with  the  estimation  of  the  percentiles.  The 

A 

simulation  study  shows  that  in  eost  cases  considered,  the  MSB  of  EDp  frost  the 
trichotoauus  response  model  is  less  than  that  of  the  dichotomous  response  model.  Although 
the  decrease  is  negligible  in  some  cases,  the  majority  of  them  are  from  5%  to  10%,  some  up 
to  20%.  Therefore,  it  is  worthwhile  to  use  the  trichotomous  response  model  if  the 
classification  of  outcomes  can  be  defined  clearly  and  implemented  easily.  Especially,  if 
one  is  interested  in  estimating  the  extreme  percentiles  and  the  test  subjects  are 
expensive,  it  would  be  very  costly  to  use  a  relatively  large  sample  sise  in  order  to  have 
reliable  estimates.  The  use  of  the  trichotomous  response  model  is  a  viable  alternative 
for  providing  better  estimates  of  these  extreme  percentiles.  Therefore  it  should  be 
seriously  considered  in  practice. 


Appendix 


Proof  of  Theorem  1 i 

Define  J  -  (B.a, . «,),  g  -  (g^gj,...^),  g'  “  (g; ,g£, . .  •  ,££>  and  let 

J(8|g)»  ate  I g*  >  be  the  Pieher  information  matrices  associated  with  the  vectors  £,£' 
respectively,  i.e.  information  from  the  (s  +  1) -response  and  s-responst  model. 

By  the  result  of  the  general  Missing  Information  Principle  (Orchard  and  Woodbury, 
1972),  we  have 


J(glg)  -  J(Slg')  +  J(g,g‘ |g>  (A) 

and  J(g>g'lg)  >  0,  which  is  the  information  lost  from  combining  0,  and  0B+1  into  one 
category.  In  particular,  J(8|fp  is  the  expectation  of  the  symmetric  matrix 


-3  iogL^i 


38 


-32iogL 


s+1 


383a. 


->2lo^1 

3B3o# 


-32logL 


s+1 


3  at 


3a* 


J 


J 


(s+1) 

j(s+1)  1 

11 

J12 

(S+1) 

j(S+1) 

21 

J22  J 

end  J(glg')  is  obtained  by  raplacing  the  last  row  and  column  of  the  above  matrix  by 


sero,  i.e.. 


J<e,R') 


Since  by  (A),  we  have  a<0 |R)  >  J(8|R'),  which  is  equivalent  to 


,(s+1)  .  , 
r11  Js 

,<s+1> 

J12 

r(s+1) 

,(s+1) 

2 1 

“22 
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which  in  turn  implies  *J22+1 '  J2*+1^  *  ®*  thus  proving  Theorem  1. 

-92logL 

The  formulae  for  are  easily  obtained  from  computing  E( — - )  snd 

- J  2 logL 

Proof  of  Theorem  2: 

Since 

-32logLj 

ae2 

-32logLj 

j(eiR')  -  e 

0 
e 

e 

0 

The  result  follows  by  applying  an  argument  similar  to  the  above. 
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Woodbury  (1972),  which  has  greatly  simplified  our  original  proofs. 
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